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1.  Introduction 


Thermoelectric  power  generation  is  the  premier  solid-state  technology  for  low- 
temperature,  greater  than  approximately  900  K,  conversion  of  heat  energy  into 
electrical  energy.  However,  the  evaluation,  interpretation,  and  analysis  of 
thermoelectric  devices  is  not  straightforward. 

The  efficiency  with  which  thermoelectric  power  generation  can  convert  heat  energy 
to  electricity  can  be  predicted,  in  part,  by  making  individual  measurements  of  the 
basic  thennoelectric  properties  of  the  constituent  materials:  the  Seebeck  coefficient, 
a;  the  electrical  resistivity,  p;  and  the  thermal  conductivity,  k.  However, 
experimental  measurement  of  those  properties  usually  results  in  significant  error, 
greater  than  ±20%.'  That  error  leads  to  significant  overestimates  of  actual  device 
performance. 

For  example,  efforts  to  quantify  the  error  of  one  property,  thermal  conductivity, 
showed  it  can  be  as  large  as  ±15%.  There  have  been  initiatives  to  characterize  these 
properties  with  significantly  improved  accuracy.2-4 

What  is  lacking  in  the  field  of  thermoelectrics  research  is  a  straightforward  method 
of  confirming  measurements  and  validating  ZT  claims.  Therefore,  one  goal  of  this 
work  is  to  demonstrate  a  simple  method  of  confirming  measurements  and  validating 
ZT  claims.  A  secondary  goal  is  the  description  of  a  straightforward  technique  for 
accurately  measuring  output  power  and  efficiency  of  thennoelectric  power 
generator  modules. 

In  this  work,  we  introduce  a  new  device  analysis  that  provides  a  simple  new 
experimental  method  to  obtain  the  usual  quantifiable  metric  for  device 
performance,  ZTmaximum.  The  significance  of  this  new  method  is  that  it  provides  a 
fast  experimental  method  to  confirm  the  validity  of  basic  materials  measurements. 
The  new  method  directly  connects  basic  materials  properties  to  conversion 
efficiency  and  employing  differential  measurements  minimizes  systemic  error. 

2.  Slope-Efficiency  Method:  Rapid  Measurement  of  Device 

ZT  maximum 


The  maximum  electrical  power  output,  Pmax,  of  any  thermoelectric  generator 
module  is  obtained  when  it  is  operated  in  an  impedance-matched  condition.  The 
impedance  matched  condition  occurs  when  internal  device  electrical  resistance, 
Rint,  exactly  equals  to  the  electrical  resistance  of  an  external  load,  Rioad.  When 
Rint  =  Rioad,  then  total  system  resistance  is  equal  to  2Rmt  and  open-circuit  voltage, 
Voc,  drops  by  half  leading  to 
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p  = 

1  mav 


4R 


int 


(i) 


For  a  thermoelectric  generator  consisting  of  some  number  “i”  of  individual 
“thennocouples”  connected  in  series  and  each  having  n-type  thermoelement  and  p- 
type  thermoelement,  the  Seebeck  effect  relates  Voc  to  the  temperature  difference, 
AT,  induced  by  the  heat  source  as  described: 

Voc  =£o(an  +  ap)iAT,  (2) 


where  an  and  ap  are  values  of  n-type  and  p-type  Seebeck  coefficients  from  each 
individual  thermoelement,  respectively.  Thus,  the  sum  of  Seebeck  coefficients  from 
i  thennocouples  is  the  ensemble-average  proportionality  between  Voc  and  AT. 
Likewise,  Rint  is  the  sum  of  resistances  from  i  thermocouples,  and  it  is  the  ensemble- 
average  electrical  resistivity  of  n-type  (pn)  and  p-type  (pp)  thermoelements  times 
their  respective  area  (A)-to-length  (f  )  values: 

Ri„.  =  a(p4+ppDf  <3) 


Thus,  P max  can  be  expressed  in  terms  of  Seebeck  coefficients: 


P  = 

A  rrmv  — 


(Xo(an+ap)i)  AT 2 
4Rint 


(4) 


This  expression  highlights  the  first  important  point:  Pmax  increases  as  AT2.  So,  for 
large  electrical  power  output,  the  largest  possible  AT  is  required. 

The  efficiency,  ®,  with  which  a  thennoelectric  generator  can  convert  heat  flow,  Q, 
to  electrical  power  is  also  important  because  the  most  electrical  power  possible 
from  a  given  amount  of  heat  flow  is  desirable.  A  new  expression  for  the  efficiency 
of  a  thennoelectric  generator  can  be  obtained  starting  with  expression  for  Pmax.  The 
ratio  of  electrical  power  generated  per  amount  of  input  heat  flow  is  the  definition 
of  efficiency: 

cp  -  (5) 


Equation  5  can  be  rewritten,  assuming  for  simplicity  a  solitary  p-n  couple  (i=l),  as 


(J>  —  vKn+Kpj  A!2 
4RjntQ 


(6) 
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The  flow  of  heat  is  dominated  by  thennal  conductivity  of  the  materials  from  which 
a  thermoelectric  generator  is  constructed,  so  Fourier’s  law  can  be  used  to  express 

Q: 


cf>  = 


(dn+otp)  AT2 
4Rint([Kn+Kp]jAT)' 


(7) 


Then  expressing  Rmt  as  described  earlier, 

^  _  _ (an+tXp)  AT2 _ 

4(PnA+PPA)([Kn+Kp]7AT) 


(8) 


For  planar  thennoelectric  generator  devices,  the  values  of  £  of  both  n-  and  p-type 
thermoelements  are  equal;  however,  cross-sectional  areas  of  n-  and  p-type  may  be 
quite  different.  Identifying  cross-sectional  area  of  n-type  as  An  and  that  of  p-type 
as  Ap  allows  a  simplification,  yielding  ®  in  terms  of  measurable  materials 
properties  and  AT : 


The  proportionality  between  ®  and  AT  is  tenned  “Zdevice” 


Z 


device 


(9) 


(10) 


Note  that  when  area-to-length  ratios  are  optimized  for  maximum  efficiency,  this 
relationship  reduces  to  the  common,  well-known  expression  for  device  ZT : 


Z 


max 


(ctn+ctp) 


C\/KnPp+ ^/KpP 


s)2 


(ii) 


The  thermoelectric  generator  efficiency  can  be  measured  as  function  of  AT,  and  the 
slope  of  that  data  should  equal  to 


d<t> 

dAT 


^device 

4 


(12) 


This  expression  highlights  a  second  important  point,  that  ®  should  linearly  increase 
as  a  function  of  AT  according  to  the  slope  indicated  by  1/4  of  the  quantity  in 
parentheses.  This  makes  sense,  because  the  ®  function  increases  linearly  with  AT, 
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and  P max  increases  as  AT2.  Taking  the  ratio  yields  a  simple  linear  dependence  on 
AT. 

It  is  important  to  note  that  the  materials  properties  are  all  temperature-dependent, 
so  taking  the  derivative  would  necessarily  yield  higher-order  tenns.  However,  we 
make  use  of  the  following  simplifications: 

1)  The  temperature  dependence  of  the  electrical  component  of  thermal 
conductivity  depends  on  the  mobility  of  charge  carriers,  and  the  electrical 
resistivity  depends  on  the  inverse  of  that  mobility,  so  these  dependencies 
cancel. 

2)  The  Seebeck  coefficient  does  have  a  relatively  small,  but  finite  temperature 
dependence;  however,  the  derivative  of  the  Seebeck  coefficient  should  yield 
a  temperature  dependence  of  T_1,  and  that  approximately  cancels  with  the 
T-1  temperature  dependence  of  lattice  thermal  conductivity  in  the 
denominator. 

The  remaining  terms  are  essentially  constant  over  the  operational  temperature 
window  of  interest  such  that  the  derivative  of  Eq.  12  would  have  a  solution  that  is 
constant.  Thus,  a  linear  slope  ought  to  be  expected.  The  linearity  of  that  slope  is  in 
fact  experimentally  observed  over  a  window  of  temperature,  as  we  show. 

A  new  index  to  determine  maximum  ZT  of  any  thennoelectric  generator  device  can 
be  obtained  by  measuring  the  slope  of  the  efficiency  as  a  function  of  AT.  To 
calculate  maximum  ZT,  4  times  the  slope  of  efficiency  multiplied  by  maximum 
temperature,  Tmaximum,  under  which  the  thermoelectric  generator  displays  linear 
behavior  with  respect  to  AT.  When  the  thennoelectric  generator  is  operated  outside 
this  linear  regime,  the  basic  materials  properties  can  no  longer  be  described  by  these 
functions  and  the  efficiency  degrades.  Therefore,  the  slope-efficiency  method 
yields  a  measure  of  maximum  ZdeviceT  as  per  Eq.  13: 

(cxn+ap)2 

(f^+ljj)(KnA„+KpAp) 

The  significance  of  this  analysis  is  it  allows  unique  means  to  rapidly  obtain 
ZTmaximum  and  confinn  properties  and  individual  measurements.  Measurements  can 
be  confinned  by  measuring  slope  of  efficiency  as  function  of  AT  and  ZT  can  be 
obtained  and  compared  to  theoretical  ZT  as  calculated  by  individual  measurements. 


-■device 


T  =  4 
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3.  Experimental  Configuration 


The  temperature  difference  across  a  thennoelectric  module,  as  defined  by  a  hot- 
side  temperature,  Thot,  and  a  cold-side  temperature,  Tcoid,  causes  electrical  power 
output  from  the  resulting  heat  flow  through  the  device.  The  magnitude  of  that  power 
is  directly  measured  by  collecting  Voc,  and  the  current-voltage  data  ranging  from 
the  Voc  condition  until  the  short-circuit  current  condition.  The  slope  of  the  current- 
voltage  curve  yields  Rint.  The  experimental  Voc  and  Rint  data  are  used  with  Eq.  1  to 
exactly  determine  the  electrical  power  output. 


To  characterize  the  thermodynamic  conversion  efficiency,  the  magnitude  of  the 
heat  flow  through  the  thennoelectric  module  must  be  known.  The  specific  strategy 
employed  in  this  work  for  determining  the  heat  flow  through  the  device  is 
schematically  shown  in  Fig.  1.  In  this  specific  strategy,  the  heat  flow  through  the 
device  is  given  as  the  sum  of  the  heat  flow  that  passes  through  a  heat  flow  meter 
after  it  leaves  the  module,  Qheat-fiow-meter,  added  to  the  magnitude  of  the  generated 
electrical  power  output,  Eq.  1 .  All  other  heat  flows  are  negligible  because  the  test 
is  performed  in  vacuum,  and  the  small  temperature  difference  between  the  cold  side 
and  the  environment  eliminates  radiative  heat  loss  at  the  cold  side. 


"hot 

Told 

Tl 


t2 


Thermoelectric 

Module 


o  Electrical 
Power 
Output 


Voc,  V(I),  and  power  measured  here 

V2 

v  nr 


maximum  electrical  power  = 


(4  Rint) 


V2 

v  n 


efficiency  = 


(4  Rint) 


V2 


Qheat-flow-meter  4  ^  ^ 


Fig.  1  Schematic  experimental  apparatus  for  measuring  the  current-voltage  data,  the 
maximum  electrical  power,  and  efficiency  of  thermoelectric  generator 


Having  eliminated  heat  flow  errors  at  the  cold  side,  the  Fourier  heat  flow  law  can 
characterize  that  heat  by  the  steady-state  temperature  difference  (T1-T2)  along  the 
length  of  a  heat  flow  meter  made  of  a  material  of  known  thermal  conductivity,  in 
this  case  oxygen-free,  high-purity  copper.  The  efficiency  is  defined,  then,  as  the 
power  output  divided  by  the  sum  of  the  Fourier  heat  flow  plus  the  power  output. 

By  increasing  the  heat  input  from  the  heat  source  in  Fig.  1,  a  range  of  temperature 
differences  across  the  thermoelectric  device  can  be  investigated.  The  temperature 
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difference  data  can  be  self-consistent  if  the  cold-side  temperature  Tcoid  is 
maintained  constant.  In  the  present  strategy,  a  constant  Tcoid  is  easily  preserved  by 
a  cooling  water  loop  at  the  base  of  the  heat  flow  meter. 

In  the  individual  cases  that  follow,  we  measure  the  power  and  efficiency  of  some 
currently  available  thennoelectric  modules.  Each  module  was  measured  using  the 
experimental  strategy  as  described.  The  obtained  data  were  then  used  to 
demonstrate  the  efficacy  of  the  new  slope-efficiency  methodology  to  analyze  and 
quantify  the  perfonnance  of  the  modules  in  tenns  of  ZT. 

3.1  Case  I:  Analysis  of  a  Precommercial  Module 

A  precommercial  module  assembled  from  unspecified  thermoelectric  materials  was 
subjected  to  the  present  analysis.  The  materials  are  designed  to  be  extremely  low 
cost,  not  high  perfonnance.  According  to  the  supplier,  the  module  technology  has 
a  hot-side  temperature  limit  of  873  K  beyond  which  slope  linearity  breaks  down, 
so  a  hot-side  temperature  of  573  K  is  sufficient  to  fully  capture  the  linear  slope  of 
interest.  In  this  test,  the  cold  side  was  maintained  at  a  constant  343  K 
(70  °C). 

Shown  in  Fig.  2  are  the  current-voltage  and  current-power  data  sets  that  were 
collected  with  Tcoid  =  343  K  with  a  maximum  temperature  difference  K  of  229  K  to 
capture  the  linear  slope.  On  the  left  panel  of  Fig.  2,  the  slopes  indicated  by  the  best- 
fit  equations  of  the  current-voltage  data  are  the  Rmt  values  as  a  function  of  the  Thot. 
The  y-intercepts  are  Voc.  The  right-side  panel  shows  the  collected  power  output 
from  the  module  ranging  from  open  circuit,  passing  through  the  parabolic  power  at 
maximum  efficiency,  and  stopping  at  short  circuit.  There  are  3  separate  AT  values 
of  135,  179,  and  229  K. 
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Current  (Amperes) 


Current  (Amperes) 

2  4 


AT(135)  =  0.6255X  - 1.9527 
R2  =  1 

AT(179)  =  0.6818X  -  2.6973 
R2  =  1 

AT(178.8)=  0.6673X  -  2.6973 
R2  =  1 

AT(228)  =  0.7046X  -  3.5732 
R2  =  1 

AT(229)=  0.7127X- 3.5827 
R2  =  1 


Fig.  2  Current-voltage  load  lines  (left)  and  current-power  curves  (right)  obtained  with 
Tcoid  =  70  °C.  The  highest  Thot  =  300  °C  for  which  AT  =  229  °C. 

The  efficiency  was  determined  and  is  plotted  in  Fig.  3.  In  this  range,  the  efficiency 
is  linear,  and  the  slope  is  0.000 1/K.  Therefore,  for  this  module,  the  determined 
ZTmaximum  0.4. 
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Fig.  3  Slope  efficiency  to  determine  ZTmax.  Tcoid  =  343  K  and  Thot  limit  =  873  K 


O  Efficiency  (%) 

-  -  Linear  (Efficiency  (%)) 

Efficiency  =  O.OOOl(AT)  -  9E-05 

R2  =  0.9996 
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3.2  Case  II:  Analysis  of  Commercial  Bismuth-Antimony- 
Selenide-Telluride  Module 

The  efficiency  measured  from  a  commercial  bismuth-antimony-selenide- 
telluride — (Bi,Sb)2(Te,Se)3 — device  is  presented  in  Fig.  4.  This  device  is  designed 
for  high  performance  for  applications  requiring  an  especially  large  value  of  thennal 
impedance,  such  as  energy  harvesting.  The  result  of  high  thermal  impedance  is  that 
the  maximum  hot-side  temperature  is  a  relatively  low  value  causing  a  relatively 
narrow  optimum  performance  window:  roughly  room  temperature  to  425  K. 

The  current-voltage  data,  current-power  data,  and  efficiency  were  collected  as  a 
function  of  AT.  The  slope  efficiency  of  0.0004/K  is  shown  in  inset  in  Fig.  4.  As 
expected,  the  slope  is  highly  linear  function  of  AT  until  deviation  from  nonlinearity 
begins  at  405  K.  The  ZTmaximum  can  be  obtained  by  the  simple  relationship,  and 
observed  maximum  temperature  of  roughly  405  K: 

ZdeviceTmaximum  —  4  (sat)  ^maximum '  (14) 

The  obtained  value  of  ZTmaximum  =  0.7  is  consistent  with  established  values  for  such 
commercial  devices  designed  for  high  perfonnance. 


Fig.  4  Slope  of  efficiency  from  (Bi,Sb)2(Te,Se)3  to  determine  ZTmax 
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3.3  Case  III:  Analysis  of  Precommercial  Lead  Telluride 
(PbTe)/TAGS  Module 

The  temperature  dependence  of  the  efficiency  of  a  precommercial  lead  telluride 
(PbTe)/TAGS  device  is  presented  in  Fig.  5.  TAGS  is  a  chemical  acronym  that 
stands  for  (GeTe)8s(AgSbTe)i5.  For  this  measurement,  the  Tcoid  was  constant  at 
300  K,  and  the  hot-side  temperature  was  increased  to  600  K.  Thus,  the  AT  across 
the  device  was  as  high  as  300  K.  For  AT  values  >150  K,  the  slope  is  small  and  a 
somewhat  nonlinear  function  of  AT.  That  temperature  range  is  outside  the  normal 
operation  window  for  PbTe/TAGS,  which  is  in  the  range  from  450  to  700  K.  A  true 
measure  of  the  representative  slope  efficiency  is  determined  where  the  efficiency 
is  linear,  greater  than  500  K. 

A  representative  value  of  ZTmaximum  can  therefore  be  obtained  by  using  Eq.  14  and 
observed  maximum  temperature  for  linear  device  behavior,  which  for  the  device 
being  measured  is  equal  to  873  K.  The  obtained  slope  is  0.0002/K  resulting  in  value 
ZTmaximum  =  0.7,  which  is  consistent  with  established  values  for  well-known 
PbTe/TAGS  modules. 


Fig.  5  Slope  of  efficiency  from  precommercial  PbTe/TAGS  to  determine  ZTma* 

4.  Conclusions 


We  have  developed  a  new,  simple  method  for  analyzing  thennoelectric  device 
performance.  The  significance  of  this  method  is  that  the  usual  quantifiable  metric 
for  device  performance,  ZTmaximum,  can  be  obtained  in  a  more  straightforward 
manner.  The  new  method  more  directly  connects  basic  materials  properties  to 
observed  conversion  efficiency  and,  employing  differential  measurements,  helps 
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minimize  systemic  error.  We  demonstrate  the  efficacy  of  this  method  by  3  separate 
cases  of  thennoelectric  power  generation  modules  fabricated  from  different 
materials.  A  low-cost  technology  demonstrates  ZT  =  0.4,  a  commercial  module 
shows  ZT  =  0.7,  and  a  precommercial  PbTe/TAGS  technology  has  ZT  =  0.7. 
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